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Dynamics of a dissipative two-level system is studied using 
quantum relaxation theory. This calculation for the first time 
goes beyond the commonly used dilute bounce gas approxi- 
mation (DBGA), even for strong damping. The new results 
obtained here deviate from the DBGA results at low temper- 
atures, however, the DBGA form is recovered at high tem- 
peratures. The results in the parameter regime 1/2 < a < 1, 
where the model has connection with the Kondo Hamiltonian, 
are of particular significance. In this regime, the spin shows a 
cross-over to a slower exponential relaxation at intermediate 
times, which is roughly half the relaxation rate at short times, 
as also observed in Quantum Monte-Carlo simulation of the 
model. The asymptotic behavior of the spin in the Kondo 
regime is in agreement with the exact conformal field theory 
results for the Kondo model. A connection of the dissipa- 
tive dynamics of the two-level system with the quantum Zeno 
effect is also presented. 

PACS indices: 67.40.Fd; 72.15.Qm; 66.30.Jt 



I. INTRODUCTION 

For the task of describing the interesting phenomenon 
of dissipation in quantum systems, two-level systems 
have served as a simple and tractable prototype. Over 
the years, several experimental situations have been 
found which can be accurately described in terms of 
damped two-level systems. Superconducting Quantum 
Interference Devices (SQUIDs) for example, are systems 
where quantum coherence can be studied at a macro- 
scopic level and description in terms of a two-level sys- 
tem is reasonably good. In metallic glasses, where atoms 
are quenched in random positions, certain atoms can end 
up finding themselves in a situation where two close by 
positions are energetically equivalent Q. So it is easiest 
for them to quantum mechanically tunnel between the 
two positions. Such a system, at very low temperatures, 
behaves like a two-level system. The sea of conduction 
electrons however impedes this motion. What one ends 
up with is a dissipative two-level system. Similar situa- 
tion exists for light interstitial particles in metals, which 
has become an interesting field of its own and goes by 
the name of "hydrogen in metals" . An example from 
the realm of quantum optics is an atom in a radiation 
field. If one is interested in the radiative properties of 
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just one excited state of the atom, one can imagine it 
to be a two-level system coupled to a dissipative "bath" 
of photons. A very well studied problem in condensed 
matter physics, that of a static spin 1/2 magnetic im- 
purity in a metal (the so-called Kondo problem) || , can 
also be thought of as a two-level system coupled to the 
conduction electrons. Later in this paper, we shall see 
some detailed application of the results to this problem. 
Recent interest in "quantum computers" , has led to the 
introduction of some interesting models to address the 
problem of quantum coherence |^^|. These too can be 
recast in the language of dissipative two-state systems. 
Lately, in the active field of high temperature supercon- 
ductivity, dissipative two-level systems have been recalled 
in the context of c-axis transport in the normal state of 
some high temperature superconductors HQ . 

A dissipative two-level system can be generally de- 
scribed by the so-called spin-boson Hamiltonian fi|. It 
consists of a psuedo-spin er, depicting the two-level sys- 
tem and a set of independent harmonic oscillators. The 
coordinates of the harmonic oscillators are linearly cou- 
pled to the z-component of the pseudo-spin a z . The 
Hamiltonian can be expressed as 

1 °° ( P 2 1 c \ 

(1) 

where A is like a tunneling matrix element if one has a 
particle in a double- well kind of a system in mind, XjS 
and pjS are the coordinates and momenta of the har- 
monic oscillators, and Cjs are their respective coupling 
constants to the two-level system, i.e., the pseudo-spin. 
Looking at the Hamiltonian in (1) one would notice that 
the (pseudo-)spin by virtue of being in a state with a z = 
+1 or -1, shifts the centers of the harmonic oscillators to 
the left or right. This results in a dissipative drag on the 
spin. It turns out that in order to analyse the influence of 
the "heat-bath" consisting of harmonic oscillators on the 
spin a, one need not demand the individual knowledge of 
CjS, rrijS and u>jS. One only requires them in a particu- 
lar combination as they appear in the so-called spectral 
density function : 

jM = ^r^^-^). (2) 

3=0 J 1 

The spectral density function contains all the information 
needed to specify the dissipative dynamics of the spin a. 
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In order to describe the dissipative behavior and intro- 
duce irreveresibilty one assumes a continuous spectrum 
of frequencies in the bath. The most common form of 
the spectral density is the "Ohmic" form where J{uS) is 
linear in u) for small frequencies: 

J{u)=auje- u}/u} % (3) 

where a is a dimensionless coupling constant parametriz- 
ing the strength of the interaction of the spin with the 
heat-bath, and lo c is a large cutoff frequency. The Ohmic 
form of the bath describes a variety of systems for various 
values of a. For example, for < a < 1/2, the spin-boson 
Hamiltonian describes the low-temperature tunneling dy- 
namics of a positively charged particle in a double-well 
potential, in the presence of conduction electrons 0] . The 
harmonic oscillators in this case represent charge density 
excitations of the electron gas, which behave like bosons. 
For the case 1/2 < a < 1, the spin-boson Hamiltonian 
mimics the Kondo model. The harmonic oscillators here 
play the role of spin density excitations of the electron 
gas, and the pseudo-spin symbolizes the magnetic impu- 
rity spin. For most purposes, the quantity of interest is 
the symmetrized spin correlation function 

c(t) = ^ < MoK(t) + <i z (ty z (o)} >, (4) 

where the angular brackets denote thermal average, and 
the Heisenberg time evolution of o~ z is dictated by the 
Hamiltonian (1). 

The spin-boson model is a very old one, but for reasons 
mentioned earlier, it still attracts considerable research 
attention. Its dynamics is what has been of most inter- 
est to researchers. It was first attacked by Leggett and 
collaborators using the path- integral formalism of Feyn- 
man and Vernon g. Using the path-integral method, it 
possible to exactly integrate over the bath degrees of free- 
dom to obtain a reduced dynamics of the pseudo-spin, at 
least formally. But in order to obtain closed expressions 
for the reduced dynamics, one has to make certain ap- 
proximations. The authors employed the so-called dilute 
bounce gas approximation (DBGA) to the functional in- 
tegral expression. The results thus obtained turned out 
to be very good for a wide range of values of the param- 
eter a and temperature. The functional integral analy- 
sis of the spin-boson Hamiltonian based on the instanton 
technique was considered by many as quite elaborate and 
formidable. This led to some simpler derivations of the 
DBGA results using second order Born approximation 
|i|.[io| and resolvent expansion techinque Jll[ . These cal- 
culations demonstrated the precise manner in which the 
DBGA is connected to conventional perturbative tech- 
niques. 

The DBGA was widely applauded for the simplicity of 
the results obtained and the accuracy with which they 
described the dynamics. For the case < a < 1/2, the 
spin shows weakly damped coherent oscillations at low 
temperatures. The physical picture is the following. The 



spin evolves quantum mechanically and is decohered by 
the effect of the heat-bath over a time which is much 
longer than A -1 . So the spin has time to evolve quan- 
tum mechanically before its coherence is destroyed. As 
the temperature is raised, the excitations of the heat-bath 
increase in strength, and at a particular temperature, the 
time over which coherence is destroyed is much smaller 
than A -1 . Here the coherent oscillations of the spin dis- 
appear completely, and what is left is an incoherent re- 
laxation. Incoherence sets in faster if a is large. When a 
becomes equal to 1/2, coherence is destroyed completely 
at all temperatures. For 1/2 < a < 1, the spin mim- 
ics a Kondo impurity spin and shows relaxation with- 
out any coherent behavior. In the region < a < 1/2, 
the DBGA is very good at not very low temperatures 
and at not very long times. In the incoherent relaxation 
behaviour, the spin decays with a rate which follows a 
power-law with temperature. At very low temperatures 
DBGA breaks down unless the coulping a is very weak. 
There have been some calculations which have gone be- 
yond the DBGA for weak coupling jl2|,[l3]]. It has been 
demonstrated that for an unbiased two-level system, for 
weak coupling, those calculations essentially yield the 
DBGA result. Thus it seems that the so-called "inter- 
bounce" interactions, at low temperature, vanish as a 2 , 
which means that for strong coupling, inter-bounce in- 
teractions have to be taken into account. This territory 
has been explored in the present investigation. In the 
Kondo regime too, the DBGA is quite good at not very 
long times, or small frequencies. In fact, DBGA yielded 
new analytical results for the dynamics of the Kondo 
spin, which reproduce many old perturbative results for 
the Kondo problem, in various limits, and generate very 
good fits to spectroscopic data in dilute magnetic alloys 
Q. Inspite of this success, DBGA fails badly at long 
times, and the question regarding the low temperature, 
long times dynamics of the Kondo spin remains open. In 
fact, behavior of C(t) in the regime T — 0, 1/2 < a < 1 
is considered a currently unresolved problem (see section 
IV E of §). 

In the present study, I put forward a calculation 
based on the resolvent expansion technique. The re- 
sults achieved, detail the dynamics of the spin at low, 
as well as high temperatures, for all relevant values of a. 
I shall present some new results in the difficult regime 
mentioned in the last para. The asymptotic behavior of 
the Kondo spin seen here agrees with the exact asymp- 
totic results derived by Affleck and Ludwig for the Kondo 
problem using conformal field theory [l5| ] . This serves as 
a stringent test of the approximation employed here. 

II. "NON-PERTURBATIVE" EXPANSION 

I begin by making a unitary transformation on the 
Hamiltonian which pulls back the centers of the shifted 
oscillators by the amount Cj(j z . The unitary operator for 
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the purpose is given by S = e 2(7 where £ = ^ c |Pi- 
The operation SHS -1 on (1) leads to 



(5) 

Now I expand the exponential factors in the first term in 
(1) in a Taylor series and write them as 



the exact treatment of Hs- This is very important from 
the point of view of the Kondo problem, because A is 
proportional to J± , the transverse part of the anisotropic 
exchange coupling in the Kondo Hamiltonian This 
modus operandi will keep both Jj_ and Ju to all orders 
and hence is expected to produce a much better pertur- 
bation theory. 

Instead of looking at the spin correlation function C(t) 
itself, it is convenient to handle its Laplace transform 
C(s): 



(±0 n 



(6) 
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Substituting (6) in (5) and regrouping terms I obtain 

t2n-l 
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Finally, I rotate my spin space basis about the y-axis, by 
7r/2 to arrive at 
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Notice that in the new basis while a x has become entan- 
gled with the bath operators, <r z is just rotated to — a x . 
Consequently, for evaluating C[t) one just has to replace 
0z by —<t Xi and 77 by H". Now the cleverest way to do 
perturbation theory would be to locate terms depending 
either only on the spin coordinates or only on the bath 
coordinates, and treat them exactly. Following this idea 
I split the Hamiltonian in (8) into a Hamiltonian for the 
spin Hs, Hamiltonian for the bath Hb, and the Hamil- 
tonian for the spin bath interaction Hj, given by 

H s = hiA<y z , 
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It is now possible to treat Hi as a perturbation, which 
does not amount to a weak coupling between the spin 
and the bath. One will notice that even if Hi is treated 
to second order, the final expression will contain a to ar- 
bitrary order due to the infinite series in £. The more 
important point is that, in addition to a, the final ex- 
pression will also contain A to arbitrary order, due to 



C(s) 



*C7(t)A, 



(10) 



where s is the Laplace transform variable. In order to 
write out the Laplace transform of (3) explicitly, the time 
evolution of o~ z in the Heisenberg representation can be 
expressed as exp(j-H"t)a z exp(—j-H"t). Introducing the 
Liouvillian associated with H as L(..) — i[H, ..], the 
Heisenberg evolution of a z can be expressed in a more 
compact notation e lLt a z . The "superoperator" e lLt com- 
pletely specifies the time evolution of the system and is 
referred to as the time evolution operator, denoted by 
U(t). Further, I assume that the initial density matrix 
of the total system is factorized as ps ■ Pb- The Laplace 
transformed correlation function then assumes the form 
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C{s) = Tracei p{a x {- 

\ 1 s — iL 



:) + (- 
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iL 



(11) 



where the (s — iL) -1 is the Laplace transform of the time 
evolution operator. As the as do not depend on the bath 
variables, one can perform a trace over the bath vari- 
ables to obtain a "bath averaged" time evolution opera- 
tor < U(s) >b , rather its Laplace transform. The tracing 
of the bath variables can be formally performed by intro- 
ducing a projection operator P which when acting on an 
operator A is defined by PA = TracesipBA). All the 
information regarding the effect of the bath on the spin is 
contained in the averaged time evolution. The major task 
then is to cvalutc < U(s) >b in a suitable approxima- 
tion. The average time evolution operator can be shown 
to satisfy the following integro-differential equation llq] 



dt 



<U(t) > B =iL s <U(t)> 1 



+ / [M(t-t)} av < U(t) > B dt, (12) 
Jo 



where [M(t — t)] av is called the memory function. The 
Laplace transform of the memory function looks like 

[M(s)] av =P(iLj){l-P) 



s-iL s -iL B -(l-P)(iL I )(l-P) 
(l-P)(iLi)P, (13) 
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where Ls, Lb and Lj are the Liouvillians associated with 
Hs, Hb and Hj, respectively. In terms of the memory 
function, the average density operator has deceptively 
compact form 



< U(a) > B - 



s - iL s - [M(s)] G 



(14) 



The form is deceptive because < U(s) >s is still exact 
and a complicated object to calculate. I shall restrict my- 
self to treating [M(s)] av to second order in L/, which is 
manageable. The corresponding expression for the mem- 
ory function takes up the subsequent form : 



[M(s)] oj; «i(PL 7 ) + (/>!,,) - 

+ P{L! 



iL s 



(PLi) 



iLs - iLs 



Li}- 



(15) 



The strategy is to first perform the trace over the bath 
states in (15), which will give a 4 x 4 matrix in the spin 
states. This matrix can then be plugged in the denom- 
inator of (14). A 4 x 4 matrix inversion then yields the 
averaged time evolution operator < U(s) >b which one 
requires. The matrix elements of the last term in (15) 
can be written in terms of correlation functions of cer- 
tain bath operators, which can be calculated using the 
properties of a bath of independent harmonic oscillators 
alone. In the ensuing analysis the spin states will be de- 
noted by Greek indices, and the number states of a set 
of harmonic oscillators by | n >, | n' > etc. We assume a 
canonical form of the density matrix at a temperature T. 
The Laplace transformed spin correlation function can 
be written in terms of the matrix elements of < U{s) >b 

C(s) = !{(+- |< U(s) > B \ +-) + (+- |< U(s) > B \ -+) 

+ (- + |< U(s) > B \ -+) + (-+ \< U(s) >b\ +-)}, 

(16) 

where | \xv) denote Liouville "states" for the spin, | ± > 
being the eigenstates of o z . One may also be interested 
in the average value of the spin <r z assuming it started 
out at t = with the value (say) o z — +1. The Laplace 
transform of < <r z (t) > is denoted in the literature by 
P(s). In the present approximation scheme, P(s) turns 
out to be identical to C(s). Another quantity of interest, 
especially for the Kondo problem, is the imaginary part 
of the spin susceptibility 



X'V) = (< + I Ps I + > - < - I Ps I - >) 

{(+- |< U(s) > B \ +-) + (+- |< U(s) >b\ -+) 

+ (-+ \< U{s) >b\ -+) + (-+ \< U(s) >b\ +-)}• 

(19) 

Some times a quantity, referred to as the spectral func- 
tion, is defined as S(u>) = x"(w)/o;. This is a symmetric 
function of u>. Now the stage is set to compute the time 
evolution operator and put the corresponding matrix el- 
ements in (16) and (19). 

III. RESULTS 

I will skip the details of the algebra involved, which 
was handled using Mathcmatica, and present the final 
form of the averaged time evolution operator. Only the 
portion of the 4x4 matrix for < U(s) >b, within the 
space | H — ), | — h) is displayed : 



1 

2^ 



-iujt 



<[a z (t),a z (0)}>dt, (17) 



which can be rewritten in terms of a Laplace transform 
as 



< U{S) >B- 



1 ( a - F{u) + F(u+) 



Det 



F(w) 



F(w) 
s - F(u) +F(lj-) 

(20) 



where 



Det = [s + F(oj+)][s + F(u-)] 

+F(u)[z+±{F(u + )+F{u>-)}]. 



(21) 



In (20) and (21) u>± = lo ±A, and F(lo) is the Laplce 
transform (with s = iu) of the following quantity, 

F(t) = A 2 exp(-2a f - e - w / wc [coth(?i/3u;) 
V Jo u 

{1 - cos(ujt)} + i sin(ujt)}doj ) + c.c. (22) 



The correlation function (16) now looks like 

C(S) = s + F{u) + [F{u + ) + F(o,_)]/2 - ^0^f 2 ' 

(23) 

The integral in (22) can be done for t, h(3 3> l/w c , so that 
its Laplace transform assumes the following structure 



F(u) 



A 2 / 2tt 



2u> c \ h/3u 



2a-l 



L(l - 2a)T(a + iujh{3/2ir) 
r(l - a + iLuh(3/2Tr) 



cos(7ra), 
(24) 



where L is Euler's Gamma function. Similarly, x"(w) 
becomes 



X"(w) - -Re 



1 r°° t tanh(^/3A/2)[ J F 1 (a; + ) - F(w_)]/2 

^J^ifejf e-<M0),a I (t)]>d t (18) [s + F( . + )][s + F( ,_ )]+fH[s+ ^ +) + F( ^ )r ^ 
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The correlation function expressed by (23) describes the 
complete dynamics of the spin in the spin-boson model. 
For T=0, one can do better than this by taking the limit 
T — > in (22) itself. In this case we do not need to make 
any approximation, and the exact form of F(u>) is given 
by 

F(u) = -i{u c /2)(A/u c ) 2 (u/Lu c ) 2a ~ 1 e^^ c T{l - 2a,Lu/u c ) 
-i(u c /2)(A/u c ) 2 (~io/uj c ) 2a ~ 1 
e -"/'- c r(l-2a > -w/w c ) > (26) 

where T (a,b) is Euler's incomplete Gamma function. In- 
terestingly co c can be completely "scaled out" of (26), 
and hence the correlation function. 



IV. DISCUSSION 

The equations (23) and (25) describe the dynamics of 
the spin at all temperatures and all values of the cou- 
pling strength. For the sake of comparison let us look at 
the corresponding DBGA result, which has the following 
form : 



C(s) = 



1 



■2F(w) 



(27) 



The reader will notice that if the A dependence in the 
arguments of the Gamma functions in (23) is ignored, 
it reduces the to the DBGA result (27). Let us discuss 
the different parameter regimes in somewhat greater de- 
tail. Figure 1 ReC(s) plotted against 10 foruj c /A = 1000, 
a = 0.1 and ksT/hA — 0.5. The solid line denotes the 
present calculation whereas the dashed line denotes the 
DBGA result. 



A. The case < a < 1/2 

For this range of values of a, the spin shows a coherent 
evolution at low temperatures. This coherent evolution is 
damped by the dissipative interaction. In this situation 
the spin-boson model can describe the coherent tunneling 
behavior of a particle in a double well potential, interact- 
ing with conduction electrons. Figures 1 shows real part 
of C(s) which is related to the structure factor for neu- 
tron scattering off the tunneling particle. The correlation 
function shows two 'inelastic' peaks which is a signature 
of coherent evolution. The position of the peaks is clearly 
shifted from the DBGA values. 

As the temperature is raised the peaks broaden and 
shift towards the origin. At a characteristic tempera- 
ture the function assumes the form of a single peak at 
lu = 0. The coherent behavior is completely destroyed 
and the spin relaxes incoherently with an exponential 



decay rate. Figure 2 shows the crossover as a func- 
tion of temperature. It is convenient to introduce a fre- 
quency A r = A(2-k /h(3uj c ) a yT(l — 2a) which is approx- 
imately the effective tunnel splitting within the DBGA. 
For T 3> A r the correlation function C(s) can be approx- 
imated by 



C(s) 



1 



s + F(0) + [F(A)+F(-A)}/2 



[F(A)-F(-A)]V4 
[F(A) + F(-A)]/2 

(28) 



The expression in (28) depicts a spin relaxing with an 
exponential decay rate which goes like T 2q_1 . The tun- 
neling partical tends to localize as temperature is in- 
creased. This is very counter-intuitive but has been ob- 
served in the several experiments on two-level systems. 
For still higher temperatures fcgT 3> fiA, one can write 
F(A) w F(— A) w -F(O). For this case the expression 
(28) reduces to the DBGA result (27). Figure 2 ReC(s) 
plotted against lu for lu c /A = 1000, a = 0.1 and various 
values of temperature r — ksT/hA. 



B. Relation to quantum Zeno effect 

One way to understand this decrease in relaxation rate 
is in terms of the so-called quantum Zeno effect. This ef- 
fect was proposed by Mishra and Sudarshan jL7| for a 
quantum system on which a series of measurements are 
made. The limit of 'continuous' measurement results in a 
freezing of the 'free' dynamics of the system completely. 
In the problem at hand, the heat-bath consisting of the 
harmonic oscillators is sensitive to the position of the 
particle, the o~ z -state of the spin in the problem at hand, 
as is obvious from the form of the Hamiltonian (1). In 
some sense the heat bath monitors the position of the 
particle in the double-well. By virtue of its posessing 
infinite degrees of freedom, it is capable of 'collapsing' 
the wave function of the tunneling particle to a particu- 
lar well. With increasing temperature the measurement 
becomes more and more efficient. In addition, the in- 
teraction being present all the time amounts to a con- 
tinuous measurement of the particle being in one of the 
two wells. Thus the coherent dynamics of the particle 
between the wells is destroyed and it tends to 'freeze' 
in one of the wells. Here one can easily work out the 
'decoherence time scale' which has been a topic of some 
recent controversy jl8 19 1, in terms of the parameters of 
the Hamiltonian. Because of the finite relaxation time of 
the particle, complete freezing will not take take place, 
and the particle will occasionally tunnel from one well to 
the other. I think this new way of looking at the dynam- 
ics of the tunneling particle helps in understanding the 
non-classical behavior seen here. 

Coming back to the issue of the dynamics of the spin- 
boson model, the value of a is crucial in deciding whether 
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the spin can have coherent dynamics or not. The coher- 
ent behavior is destroyed at a = 1/2. At a = 1/2, C(t) 
decays with a single relaxation rate ttA 2 /2lo c , which is 
what DBGA also yields. 



C. The case 1/2 < a < 1 

Beyond a = 1/2 we enter the interesting regime where 
the spin-boson model is related to the Kondo Hamilto- 
nian. As mentioned earlier, the results (23)- (25) contain 
A, and hence J±, to arbitrary order. Hence, the calcula- 
tion is expected to be better than that in [[14j . In order to 
obtain exponential relaxation rates which are relevant for 
experimentalist, I first neglect the frequency dependence 
in the arguments of the Gamma functions in (23). The 
function C(s) is anyway concentrated around lo =0. This 
results in the following form of the correlation function 



C(s) 



s + F(0) + [F(A)+F(- 



A\l in , [F(A)-F(-A)]V4 
-^)l/ Z + s+[F(A)+F(-A)}/2 

(29) 



The above expression can be decomposed into partial 
fractions as 



C(s) 



D 



s + r 1 s + r 2 



(30) 



where 



_[F(0)+F(A) + F(-A)] 
±^ 2 (0)-[F(A)-F(-A)]2, 



A,B = -± 



F(0)/2 



2 y /F 2 (0)-[F(A)-F(-A)] 2 



(31) 



The spin shows two exponential decay rates Ti and IV 
The Quantum Monte Carlo simulation of the spin-boson 
model by Egger and Weiss revealed that after an initial 
exponential relaxation, the spin shows a slower exponen- 
tial relaxation p0| ]. This second slower decay rate was 
estimated to be roughly half the initial decay rate. The 
initial decay rate is seen to match the DBGA value. In 
order to make contact with the QMC prediction I con- 
sider the special case F 2 (0) > — [F(A) -F(-A)] 2 which 
is satisfied for the parameters chosen by Egger and Weiss. 
For this case I arrive at 

ri « F(0) + [F(A) + F (-A)]/2 
T 2 « [F(A) + F(-A)]/2 

-[F(A)-F(-A)] 2 /2 



B 



F 2 {0) - [F(A) - F(-A)] 2 
T (A)-F(-A)] 2 /2 



F 2 (0) - [F(A) - F(-A)} 2 



(32) 



Clearly in (32), T 2 ~ Ti /2. Moreover the quantity B = 
— [-F(A) — F(— A)] 2 being small, Tx dominates for small 
times. The second decay being slower, shows up at long 
times. 

In the Kondo regime, expression (23) for C(s) is valid 
for all temperatures of practical interest. Thus the dy- 
namics is also described at very low temperature and 
small frequency, a regime inaccessible in 14 1. To stress 
this point let us look at the absorptive part of the spin 
susceptiblity x"(lS). An exact relation derived by Shiba 
pl[ for the general Anderson model says that at T = 0, 
x(lo) should vanish linearly in lo : 



(33) 



This implies that C(t) will asympotically behave as ~ t~ 2 
||. DBGA predicts an asymptotic behavior ~ t^ 2 ^^ a \ 
and hence a divergent spectral function S(lu). Very re- 
cently, conformal field theory has been employed to calcu- 
late the exact asymptotic spin correlation function for the 
Kondo Hamiltonian |t5) . The exact results show that the 
spin correlation function asymptoticaly decays as ~ t~ 2 . 
Let us take a closer look at the x"(w) calculated here. 
In order to make the complicated expression (25) more 
tractable, I consider the case u> < ksT/Ti <C A. In this 
situation F(u±) become temperature independent, and 
F(u)) becomes negligible in their comparison. It is conve- 
nient to introduce a parameter e = 1 — a which is related 
to the Jii in the Kondo problem. In reality, it may be 
treated as a small parameter. For small lo one may ex- 
pand the various terms involved in a Taylor series : 



F(u+) +F(w_) 
F(u>+) - F(u-) 

where 
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2 7 A 1_2e sin(Tre) + 2i7wA _2e cos(?re) 
2i~fA 1 ~ 2e cos(Tre) + 2jljA^ 2c sin(Tre) (34) 



(35) 



4u c r(2e) sin(Tre) ' 
In this approximation x"( w ) assumes the following form 
1 



-Re 



iA + lo tan(7re) 



r iujA 2fL 
L 7 cos(7re) 



A tan(7re) 



A 2 



7 cos(7re) [iA + lo tan(7re)] 2 
(36) 



From the above expression it is clear that, 
lim w _»o ~ lo. This implies a finite static suscep- 

tiblity, and an asymptotic correlation function going as 
~ t~ 2 , which is a rigorous result from conformal field 
theory p!q |. The major achievement of this work is that 
the low temperature and long time dynamics of the spin 
in the Kondo regime is solved and the results seem quite 
satisfactory. 

In conclusion, I have presented a calculation of the 
dynamics of a strongly damped two-level system which 



G 



goes beyond the DBG A. In the coherent regime the re- 
sults substantially deviate from the corresponding DBGA 
ones. In the Kondo regime, the short and intermedi- 
ate time behavior agrees with the Quantum Monte Carlo 
simulation of Egger and Weiss. The spin crosses over to 
a slower exponential relaxation after an initial fast decay. 
The asymptotic behavior of the spin is in agreement with 
the exact conformal field theory results. An interesting 
relation to the quantum Zeno effect is presented. Finally, 
it should be mentioned that it is straightforward to carry 
out a similar calculation for a biased two-level system. 
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